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ABSTRACT 

Magnetic (or tidal) interactions between "hot Jupiters" and their host stars can potentially enhance 
chromospheric and coronal activity. An ideal testbed for investigating this effect is provided by 
the extreme WASP-18 system, which features a massive (~10 times Jupiter) close-in ( < 1 day period) 
transiting planet orbiting a young F6 star. Optical and X-ray observations of WASP-18 were conducted 
in November 2011. The high-resolution echelle spectrograph MIKE was used on the 6.5m Magellan 
Clay Telescope to obtain 13 spectra spanning planetary orbital phases of 0.7-1.4, while the X-ray 
Telescope on Swift provided contemporaneous monitoring with a stacked exposure of ~50 ks. The 
cores of the Ca II H and K lines do not show significant variability over multiple orbits spanning ^8 d, 
in contrast to the expectation of phase-dependent chromospheric activity enhancements for efficient 
star-planet interaction. The star is also X-ray faint, with logLx < 27.6 erg s _1 (0.3-2 keV), indicating 
that coronal activity is likewise low. The lack of detectable star-planet interaction in this extreme 
system requires that any such effect must here be transient, if indeed present. We demonstrate that 
searches for Ca II H and K variability can potentially mistake a stellar hotspot, if observed over a 
short segment of the rotation period, for planet-induced activity. Taken together, these results suggest 
that the utility of star-planet interaction as a robust method of estimating exoplanet magnetic field 
strengths may be limited. 

Subject headings: planetary systems — stars: activity — stars: individual (WASP-18) 



1. INTRODUCTION 

The possibility of magnetic star-planet interactions 
was initially explored by Cuntz et al. (2000; see also 
Rubenstein & Schaefer 2000), and this and subsequent 
work has indicated that reconnection events could be an 
important effect in "hot Jupiter" systems, acting to pro- 
duce enhancement^ in chromospheric and coronal ac- 
tivity co-rotating with the planet rather than with the 
stellar rotation. It has been conjectured that the en- 
ergy available from reconnection events should scale with 
the product of the stellar and planetary magnetic field 
strengths, as well as the relative velocity of the mag- 
netic field lines, and inversely with approximately the 
square of the orbital semi-major axis (Cuntz et al. 2000; 
Kashyap et al. 2008). Consequently, a calibrated rela- 
tionship between the amplitude of star-planet interaction 
and stellar and orbital parameters could permit estima- 
tion of exoplanet magnetic field strengths. Throughout 
this work we quantify interaction strength in otherwise 
similar systems as Mp/a 2 , where a is the semi-major 
axis in AU and the planetary magnetic field strength is 
taken to scale with the planetary mass Mp (e.g., Arge et 
al. 1995; Stevens 2005); this is broadly consistent with 
the various proposed trends investigated (e.g., Kashyap 
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et al. 2008; Shkolnik et al. 2008; Scharf 2010; Poppen- 
haeger et al. 2010). Numerical studies suggest that star- 
planet interactions can potentially generate sufficient en- 
ergy to be observable as, e.g., phase- variable core Ca II 
H and K or X-ray emission (probing chromospheric and 
coronal activity, respectively) in monitoring of individual 
hot Jupiter systems, or as a greater average level of ac- 
tivity in systems with more massive or closer-in planets. 
For example, Lanza et al. (2008, 2011) modeled chro- 
mospheric hot spots in several systems (offset from the 
subplanetary point by varying degrees) as arising from 
star-planet magnetic reconnection events, and Cohen et 
al. (2009, 2011) carried out three-dimensional magneto- 
hydrodynamic simulations demonstrating that close-in 
giant planets can produce an increase in overall stellar 
activity and generate (non-persistent) coronal hot spots 
that rotate synchronously with the planet (albeit poten- 
tially shifted in phase); see also Pillitteri et al. (2010). 

Observational evidence of magnetic star-planet inter- 
action has now been presented for several individual sys- 
tems. For example, Shkolnik et al. (2005) examined 10 
stars (Kl to F7) known at the time to possess mas- 
sive, close-in planets (median minimum planetary mass 
Aipsini — 0.6Mj up , median orbital period P = 3.4 d), 
in low-eccentricity orbits, and claimed evidence of star- 
planet interaction in HD 179949 and v And (both F8) 
based on slight Ca II H and K emission variability syn- 
chronized with the planetary period (see also discussion 
in Gu et al. 2005), although this synchronization is ap- 
parently transient (Shkolnik et al. 2008). Chandra ob- 
servations of HD 179949 showed variable X-ray emission 
with an apparent maximum near the phase associated 

6 The cores of the Ca II H and K lines are good indicators of 
chromospheric activity; residuals from Ho or the calcium infrared 
triplet are also useful for this purpose. 
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with Ca II H and K enhancement (Saar et al. 2008). 
Fares et al. (2012) carried out spectropolarimetric obser- 
vations of HD 179949, finding that the stellar magneto- 
sphere is highly tilted (producing two maxima per rota- 
tion period) and that chromospheric activity is primarily 
linked to stellar rotation, although low-level fluctuations 
near the beat period could be planet-induced (see also 
Gurdemir et al. 2012). Later Chandra observations of 
v And (as well as concurrent optical spectroscopy) did 
not find any indications of star-planet interaction (Pop- 
penhaeger et al. 2011a). Interaction was claimed for 
r Boo (F7) based primarily on photometric spot mod- 
eling, although the tidal-locking of the star to the planet 
complicates interpretation (Walker et al. 2008; see also 
Shkolnik et al. 2008). HD 189733 (K0) twice showed an 
X-ray flare after secondary eclipse (Pillitteri et al. 2010, 
2011) but chromospheric variability appears to be tied to 
stellar rotation (Shkolnik et al. 2008; Fares et al. 2010). 
Optical spectroscopy of several systems (including v And 
and t Boo) carried out by Lenz et al. (2011) did not find 
evidence of planet-induced chromospheric activity, but 
did identify potential interaction between the M dwarf 
HD 41004B (which has binary, planet-bearing K dwarf 
companion HD 41004A) and its brown-dwarf companion 
HD 41004Bb. As a counterexample, 51 Peg (G5) is one 
of several systems in which a sensitive search did not un- 
cover evidence of star-planet interaction; its low Lx and 
low logi? HK suggest it may in fact be in an extended 
state of especially low activity, analogous to the Sun's 
Maunder minimum (Poppenhaeger et al. 2009). 

Several other studies have used large samples of planet- 
bearing stars to test for significantly increased activity 
in hot Jupiter systems, statistically averaging over or- 
bital phase. X-ray investigations of potential planet- 
induced enhancements in coronal activity have to date 
provided mixed results. This may be due to the diffi- 
culty of controlling for selection effects (more active stars, 
which tend to be brighter in X-rays, are disproportion- 
ately identified as planet-bearing with high mass, close-in 
planets, since those stellar properties limit detectability 
of low mass, distant planets) or may be related to sam- 
pling a heterogeneous mix of late- type (typically F, G, 
K, and M) stellar classes. Kashyap et al. (2008) con- 
sidered a large sample primarily drawn from ROSAT 
pointed and all-sky survey (RASS) observations (con- 
taining ^30%/70% X-ray detections/upper limits) and 
found that, after attempting to control for selection ef- 
fects, stars with close-in planets (a < 0.15 AU) were 
1.3-4 times more X-ray luminous than stars with dis- 
tant (a > 1.5 AU) planets. Similarly, Scharf (2010), us- 
ing RASS data again with a high percentage of upper 
limits, did not find enhanced X-ray emission for close- 
in planets (note that their distant planets encompassed 
a > 0.15 AU) and estimated their Lx — Mp correlation 
to be robust against observational bias. On the other 
hand, Poppenhaeger et al. (2010) looked at a volume- 
limited (d < 30 pc) sample of stars observed with XMM- 
Newton or ROSAT and found no significant correlations 
of Lx/Lb i with Mp or a; they explain the possible cor- 
relation of Lx with Mp/a as arising entirely from se- 
lection effects. Importantly, Poppenhaeger et al. (2011b) 
demonstrated that distance in shallow RASS datasets co- 
varies with both Lx and Mp (indicating that deeper ob- 
servations are essential to sidestep selection effects) and 



that apparent correlations with Lx vanish when Lx/L\, Q \ 
is used instead (thereby controlling, to first order, for 
spectral class) . Large-sample studies of potential planet- 
induced enhancements in chromospheric activity have 
also given mixed results; Canto Martins et al. (2011) 
find no obvious correlations between planetary mass or 
semi-major axis and the Ca II H and K activity indicator 
log\R HK , whereas Krejcova & Budaj (2012) do find these 
variables to correlate but only within the subset of cooler 
host stars with T e g < 5500 K. 

In this paper, we present a sensitive search for star- 
planet interaction in the extreme WASP-18 system. If 
star-planet interaction is to become a useful probe of 
exoplanet magnetic field strengths in hot Jupiter sys- 
tems, additional convincing instances, beyond the hand- 
ful suggested to date, must be identified. Observa- 
tions of a strongly-interacting system could constitute 
a contextual template to guide interpretation of results 
in more weakly-interacting systems (including those al- 
ready studied). Additionally, a robust measurement of 
star-planet interaction in an extreme high-mass, short- 
period system would supply productive leverage for un- 
covering and later quantifying scaling relations. Con- 
versely, a failure to detect planet-induced stellar activity 
in an extreme system would severely constrain the practi- 
cal relevance of star-planet interaction and could impact 
theoretical understanding and numerical modeling of this 
effect. 

1.1. WASP-18 in context 

We identified WASP-18 as the potentially most 
strongly interacting system among currently discovered 
exoplanets, a distinction it retains as of April 2012. We 
obtained a list of all confirmed planets with Mp > O.lMj 
(i.e., planetary mass greater than 10% that of Jupiter), 
orbital period P < 100 d, and V < 15 from the Exo- 
planet Orbit DatabaseQ (Wright et al. 2011). As may be 
seen in Figure 1, WASP-18b (Hellier et al. 2009; South- 
worth et al. 2009), with M P = 10.4Mj, P = 0.941 d, 
and a — 0.020 AU, has the largest value of the inter- 
action strength proxy Mp/a 2 within this sample (and 
would also rank first for other plausible parameteri- 
zations of interaction strength). More specifically, af- 
ter excluding brown dwarfs (with Mp > 13Mj up ), of 
the 551 entries in the Exoplanet Orbit Database, the 
Mp/a 2 ~ 25000 M Jup AU -2 for WASP-18 is nearly 2.5 
times greater than the next largest value (CoRoT-14b), 
and only 5 systems have Mp/a 2 > 5000 Af Jup AU~ . 
Because WASP-18 is a transiting system, the plane- 
tary properties are securely established by lightcurve and 
radial- velocity measurements. The values of Mp/a 2 are 
also labeled for several systems described in §1; for ex- 
ample, HD 179949 has Mp/a 2 ~ 500 M Jup AU" 2 . 

The stellar and planetary properties for the WASP-18 
system are given in Table 1, as are those for previously 
identified star-planet interaction candidates HD 179949, 
v And, and r Boo (but recall r Boo is tidally locked, 
with stellar and orbital periods of ^3.2 d, which may 
suppress magnetic interaction). The value of Mp/a 2 for 
WASP-18b exceeds that for the other candidate systems 
by 1-2 orders of magnitude. Such extreme planetary 
systems are difficult to maintain, and WASP-18b itself is 

7 Available at http://www.exoplanets.org 
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Figure 1. Distribution of star-planet interaction strength, assumed here to scale as Mp /a 2 , plotted versus apparent V magnitude 
(top) and star-planet separation in units of stellar radii (bottom). Note that the horizontal axis is logarithmic. Transiting planets 
are plotted in blue, and symbol size scales as Mp' 6 . Brown dwarfs (> 13Mj) are plotted as open symbols. Radial-velocity masses 
are lower limits, but establishing sin i would typically result in only modest rightward movement of those points. Several relevant 
or interesting systems are labeled (the b planet designation is implicit). 



expected to have a short lifetime against infall from tidal 
drag (Hellier et al. 2009). The star WASP-18 has spec- 
tral type F6, effective temperature T e g — 6400 ± 100 K, 
and stellar mass M* = 1.22 ± 0.03 M© (values from 
the Exoplanet Orbit Database, from which original ref- 
erences may be obtained) ; conveniently, these are similar 
to the properties of HD 179949 (F8, 6170 K, 1.18 M ), 
v And (F8, 6210 K, 1.31 M Q ), and r Boo (F7, 6390 K, 
1.34 M©). WASP-18 is slightly hotter than HD 179949 
and v And, and consequently likely possesses a slightly 
shallower outer convection zone. On the other hand, 
the stellar rotation period is shorter, at ^5.6 d, esti- 
mated from usini = 11 km s _1 ; the orbit and stellar 
rotation are aligned^ (Triaud et al. 2010). HD 179949 
and v And have rotational periods of 7.6 d (Fares et 

8 All four of these systems are near the T e fj ~ 6250 K border 
above which hot-Jupiter hosts tend to have high stellar obliquities 
(Winn et al. 2010), but the spin-orbit alignment in WASP-18 is 
more characteristic of cooler hot- Jupiter hosts. 



al. 2012) and -12 d (Shkolnik et al. 2008), respectively. 
Further, WASP-18 appears to be a younger star, with 
an age of ~500-700 Myr (Hellier et al. 2009; Brown et 
al. 2011). The age-rotation-activity relation would then 
predict relatively greater intrinsic activity in WASP-18, 
apart from any planet-induced modulation or enhance- 
ment (but see §3.2 for caveats). Observationally, WASP- 
18 has other appealing features: the star is among the 
brighter transit-detected systems, and the short plane- 
tary period facilitates rapid accumulation of phase cov- 
erage across multiple orbits. Due to its extreme proper- 
ties, WASP-18 has been studied at a range of frequencies; 
for example, Nymeyer et al. (2011) used Spitzer obser- 
vations of secondary eclipse to infer that the planet has 
T ~3100 K with near-zero values for both albedo and 
day/night side energy redistribution. However, high- 
resolution, high signal-to-noise Ca II H and K spec- 
troscopy and sensitive X-ray observations have not been 
published prior to the observations presented here. 
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2. OBSERVATIONS 
2.1. Data acquisition and reduction 

Optical spectroscopy was carried out with the 6.5m 
Clay Telescope on Nov. 6-7 and 10-13 2011 using the 
Magellan Inamori Kyocera Echelle (MIKE), a high- 
throughput double cchcllc spectrograph. A 0.7x5" slit 
was used throughout, providing a resolution of R~40000 
near the Ca II H and K lines. A communication board 
failure prevented use of the blue-side CCD on Nov. 10 
and bad weather limited observing on Nov. 11. In to- 
tal, 13 visits to WASP-18 were obtained, with combined 
exposure times of 20-40 minutes per visit. The plane- 
tary orbital phase coverage (Figure 2) spans 0.7-1.4 and 
includes one visit at central planetary transit. 

The spectra were reduced using the IDL MIKE Redujfl 
package. Briefly, milky flats (taken of 5 Ori with the dif- 
fuser) were combined and normalized to correct pixel re- 
sponse variations; trace flats (from internal lamps) were 
combined for order and slit tracing; arc images (includ- 
ing one associated with each science exposure) were pro- 
cessed and used to derive two-dimensional wavelength 
maps, using a fifth-degree polynomial fit along the orders 
which was found to give generally structureless residuals; 
the slit profile was determined for each order; science ex- 
posures were processed (bias-subtracted from overscan 
and flat-corrected), the object was traced and the sky 
subtracted, and the object flux was optimally extracted. 
Flux is retained in relative units (i.e., no spectrophoto- 
metric standard star calibration was conducted) as we 
are interested in normalized spectra. Automated cosmic 
ray correction was not performed due to concern about 
potential introduction of spurious variability; two science 
frames with an obvious cosmic ray located near the two- 
dimensional center of the Ca II K line were discarded. 

Contemporaneous monitoring was obtained with the 
Swift satellite from Nov. 1-17 (Target ID 32149), no- 
tably with the X-ray Telescope (XRT) through 58 snap- 
shots of 500-1500 seconds each. The XRT observations 
were reduced using HEASARC version 6.11.1. An ex- 
posure map was generated for each sequence, including 
a vignetting correction and with time-dependent instru- 
ment maps used whenever attitude variations exceeded 
2.4". Images were extracted in both the full (0.3-10 keV) 
and soft (0.3-2 keV) bands. The stacked XRT effective 
exposure from all sequences at the position of WASP- 
18 reached ~50 ks. The UV/Optical Telescope (UVOT) 
was also used to obtain both photometry (with the U, 
UVW1, UVM2, and UVW2 filters) and UV grism spec- 
troscopy, with coverage of the Mg II h and k lines which 
are also chromospheric activity indicators (Buccino & 
Mauas 2008). 

2.2. Analysis and results 

We checked for phase-dependent variability, as from 
one-sided enhanced chromospheric activity, within the 
cores of the Ca II H and K lines, focusing primarily 
upon the latter. The Ca II K line and nearby stellar 
continuum were isolated through selecting wavelengths 
spanning 3920-3950A, which included contributions from 
two orders. For each visit, the flux from individual ex- 
posures was summed, and the wavelength scale was ad- 

9 http : //web .mit . edu/~bu rles/www/MIKE/ 
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Figure 2. Phase coverage of the Magellan observations, plot- 
ted on the stellar radial velocity curve from the ephemeris of 
Hellier et al. (2009; error here < 0.001). Dates are UT 2011; 
phase coverage of the 0.94 d planetary orbit was accumulated 
over five nights. The Nov 12 and first Nov 14 points have 
slight vertical offsets for clarity. The vertical dotted lines 
mark planetary transit. 

justed through subtraction of a center-of-mass velocity of 
3.1961 km s _1 and a stellar reflex velocity with an am- 
plitude of 1.8183 km s _1 , using the values and ephemeris 
from Hellier et al. (2009). Continuum normalization near 
the Ca II H and K regions is difficult because the absorp- 
tion lines are quite broad relative to the blaze function, 
there are a large number of nearby metal lines, and the 
Wien tail is falling off steeply. We used a similar normal- 
ization method to that conducted and detailed by Shkol- 
nik et al. (2005, 2008), removing a slight linear trend 
over ~7 A centered near the line cores, such that unity 
represents approximately one-third of the local stellar 
continuum. The normalized Ca II H and K spectra for 
all 13 visits are plotted in Figure 3 (top). 

Residuals from median normalized Ca II H and K 
spectra were calculated for each visit, and then stacked 
(weighted by exposure times corrected for seeing) within 
orbital phase bins of 0.70-0.93, 0.99-1.19, and 1.21-1.41; 
these bins then contain spectra from 4, 5, and 4 visits, 
respectively. The residuals were boxcar smoothed with a 
width of 0.4A and then any remaining slight linear trend 
was subtracted to avoid influencing comparison of the 
line core. There is no obvious difference in the resid- 
ual spectra within the three separate phase bins, to an 
amplitude of 0.003 of the normalized continuum (Fig- 
ure 3, bottom). Integrating the normalized flux over the 
1 A line cores within phase bins of 0.70-0.93, 0.99-1.19, 
and 1.21—1.41 gives nearly identical values: 0.252±0.001, 
0.251±0.002, and 0.253±0.003 for K and 0.225±0.001, 
0.227±0.001, and 0.228±0.003 for H, where the la un- 
certainties have been estimated as the standard deviation 
of the component flux measurements within each phase 
bin. The 3a limit on the relative flux variability is <3.2% 
and <3.4% for the K and H line cores, respectively. The 
apparent marginal tendency for the 1.21-1.41 K bin to 
show relatively more red-side emission is within the noise 
leveO but in any case this slight profile variation is not 
present in H, is lopsided rather than symmetric about the 

10 This may be verified through comparison to the residuals 
near photospheric absorption features, which should not show any 
variability phased with the planetary orbit. 
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Figure 3. Normalized spectra of the Ca II H (right, top) and K (left, top) regions overplotted for each of 13 visits. There is no 
obvious change in the depth of the central core of the Ca II K line. Residuals from the median spectrum, binned by planetary 
orbital phase and boxcar smoothed with a width of 0.4 A, are shown at bottom (after removal of slight linear trends). The 
vertical dotted lines mark the Ca II H and K cores. 



core, and lags rather than leads the subplanetary point; 
these effects are inconsistent with planetary magnetic in- 
teraction. 

The XRT observations were initially planned for simi- 
lar phase-resolved analysis, but the star proved to be un- 
expectedly X-ray faint (see also §3.2) and so all sequences 
were combined for an effective exposure of ^50 ks. 
WASP- 18 is not detected by XRT in the stacked ex- 
posure, in either the soft (Figure 4) or the full bands. 
Specifically, within an aperture of 20" centered on the 
SIMBAD optical position of WASP-18, there are 3.8 
counts in the 0.3-2 keV band, where the expected back- 
ground is 3.7 counts. The 95% upper limit on the source 
counts within 20" is < 5.9 net counts (from Kraft et 



al. 1991), corresponding to < 8.1 total net counts after 
accounting for the XRT point spread function (Moretti 
et al. 2005). For a coronal model with kT=l keV and 
solar abundances the unabsorbed 0.3-2 keV X-ray flux 
is < 3.3xl0~ 15 erg s _1 , calculated using the Portable 
Interactive Multi-Mission Simulator (PIMMSE3) for a 
plasma/APEC model presuming an intervening column 
of -/Vh = 10 18 ~ 19 cm -2 ; other plausible models give simi- 
lar results to within < 30% (or 0.1 dex). At the Hipparcos 
distance to WASP-18 of 100 ± 10.6 pc, this corresponds 
to a limiting X-ray luminosity of logLx < 27.6 erg s _1 
(0.3-2 keV). 



http : // cxc . harvard . edu/ toolkit /pimms . j sp 
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3. DISCUSSION 

We compare the observed persistent low activity in 
WASP-18 to expectations for planet-induced variability, 
as well as to the intrinsic properties of similar stars, 
and additionally explore the general possibility of stel- 
lar hotspots acting to mimic star-planet interaction. 

3.1. Expected variability in WASP-18 

Past work would seem to suggest that the degree of 
planet-induced Ca II H and K variability expected for 
WASP-18 could be substantial. For example, Shkolnik et 
al. (2005, 2008) report "on/off" Ca II H and K variability 
phased with the planetary orbital period in HD 179949 
and v And (see also Poppenhaeger et al. 2011a), which 
have stellar types of F8 and do not show obvious core 
emission within the deep photospheric absorption (sim- 
ilar to WASP-18). The amplitude of residual flux asso- 
ciated by Shkolnik et al. (2008, Figure 3; 2005, Figure 
6) with star-planet activity in HD 179949 (v And) is 
-0.017 (-0.008)0 The value of M P /a 2 for WASP-18 
could suggest an effect 50 (130) times greater than for 
HD 179949 (v And) 03 In contrast to such predictions, 
the level of phase-binned variability observed for WASP- 
18 is < 0.003 (Figure 3; §2.2), below the levels detected 
in HD 179949 and v And. 

The degree of X-ray variability expected in WASP- 
18 is more difficult to estimate, as only a handful of 
X-ray brightening events have been attributed to po- 
tential star-planet interaction. For example, HD 179949 
displayed an increase in X-ray emission by ^30% near 
the phase associated with the Ca II H and K enhance- 
ment (Saar et al. 2008); with a stellar X-ray luminos- 
ity of logLx = 28.6 erg s" 1 (Kashyap et al. 2008), this 
is an increase of logLx = 28.1 erg s _1 . HD 189733 
showed two X-ray flares near phase 0.53 (i.e., shortly 
after occultation) , peaking at twice the baseline count 
rate (Pillitteri et al. 2011); with a stellar X-ray lumi- 
nosity of logix = 28.4 erg s -1 (Kashyap et al. 2008), 
this is an increase of logLx = 28.4 erg s -1 . Cohen et 
al. (2011) used their MHD modeling of the HD 189733 
system to estimate that the energy available from mag- 
netic reconnection to accelerate particles into the stellar 
corona is ~10 28 erg s _1 (after applying conservative ef- 
ficiency assumptions). From these examples, WASP-18 
might be expected, taking into account its exceptional 
value of Mp/a 2 , to possess interaction energy available 
for coronal activity enhancement sufficient to generate 
absolute increases in its X-ray luminosity by several times 
10 28 erg s _1 . However, the observed XRT upper limit of 
logLx < 27.6 erg s _1 indicates that the X-ray luminosity 
did not achieve this level for any more than < 10% of the 
Swift coverage. Without even a stacked X-ray detection 
it is not possible to evaluate the relative amplitude of 
phase-dependent X-ray variability, if any, in WASP-18. 

The scaling of star-planet interaction strength with 

12 The ratio of residuals for HD 179949 and v And is ~2, while 
the ratio of Mp/a 2 is ~2.5. However, given the dissimilarity of 
other parameters, such as intrinsic stellar activity, and the apparent 
variability of any star-planet interaction in these systems, we do 
not consider this necessarily informative. 

13 While less physically motivated, alternative scalings with ei- 
ther separately 1/a 2 or Mp would still predict a substantially 
greater effect for WASP-18, by a factor of 4.8 (8.7) or 11 (15) 
compared to HD 179949 (v And). 
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Figure 4. Top: Digital Sky Survey image of WASP 18. Mid- 
dle: Swift XRT 0.3-2 keV image of WASP 18 region, Gaus- 
sian smoothed with a kernel of 7 pixels and shown with log- 
arithmic scaling. The green circle has a radius of 20" and 
is centered on the SIMBAD position of WASP 18. Bottom: 
XRT exposure map; the color bar gives the effective exposure 
in seconds. 

Mp/a 2 that we (and others; §1) have adopted approx- 
imately holds for otherwise similar systems. We briefly 
consider the impact upon the energy available for star- 
planet interaction due to potentially differing planetary 
or stellar magnetic field strengths in WASP-18 com- 
pared to HD 179949 and v And. The planetary mag- 
netic field strength depends most relevantly upon Mp 
but also scales inversely with the planetary rotational 
period. For tidally-locked planets, as all these are be- 
lieved to be (out to ^0.15 AU; Bodenheimer et al. 2001), 
the planetary rotation period is identical to the orbital 
period, and both are substantially shorter for WASP-18 
than for HD 179949 or v And. This suggests (Arge et 
al. 1995; Sanchez-Lavega 2004; Stevens 2005) an increase 
in baseline planetary magnetic field strength by a factor 
of a few, which would enhance any star-planet interac- 
tion. The relative velocity between magnetic field lines, 
which is governed by the difference between the plane- 
tary and stellar rotation rates, is v rc \ = K(R*/a)—v I0 t for 
tidally-locked planets, where K is the orbital velocity and 
w ro t is the equatorial stellar rotational velocity (Cuntz et 
al. 2000). This would suggest a further increase in inter- 
action strength by a factor of a few. On the other hand, 
the interaction strength approximately scales with the 
stellar magnetic field, B t (Cuntz et al. 2000; Kashyap et 
al. 2008; see Lanza 2009 and Scharf 2010 for slightly dif- 
ferent formulations), which appears to be substantially 
weaker in WASP-18 than in at least HD 179949 and 
perhaps also v And. This is evident by the X-ray non- 
detection of WASP-18, as is observed to depend ap- 
proximately linearly upon Lx (Pevtsov et al. 2003). That 
log (Lx/^boi) is at least one order of magnitude lower 
in WASP-18 than HD 179949 possibly offsets any gains 
in interaction strength due to faster planetary rotation, 
although the difference in inferred stellar magnetic field 
strengths is less between WASP-18 and v And. However, 
since the observed amplitude of Ca II H and K variabil- 
ity in WASP-18 is already > 5 times lower than that sug- 
gested to arise from star-planet interaction in HD 179949, 
B* would need to be > 50 times lower in WASP-18 than 
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Table 1 

Candidates for star-planet interaction 









Star 






Planet 








Activity 






Name 


Type 


V 


B-V 


d 


M P 


P 





Age 


Mp/a 2 


^HK 


logLx 












(PC) 


(Mj) 


(d) 


(AU) 


(Gyr) 


(Mj/AU 2 ) 




(erg s" 1 ) 




WASP-18 


F6 


9.39 


0.49 


100 


10.4 


0.94 


0.0201 


0.6 


25700 


<-5.0 a 


<27.6 


<-6.2 


r Boo b 


F7 


4.50 


0.51 


15.62 


4.1 


3.31 


0.0480 


2.5 


1780 


-4.73 


28.8 


-5.3 


HD 179949 


F8 


6.25 


0.55 


27.5 


0.95 


3.09 


0.0439 


2.1 


490 


-4.80 


28.6 


-5.3 


v And 


F8 


4.10 


0.54 


13.49 


0.69 


4.62 


0.0594 


3.8 


200 


-5.07 


27.6 


-6.5 



Note. — Data is taken from the Exoplanet Orbit Database (http://www.exoplanets.org), except for planetary mass and system age 
which are from the Extrasolar Planets Encyclopedia (http://www.exoplanet.eu), and logLx which is from this work, Poppenhaeger et 



al. (2012), Kashyap et al. (2008), and Poppenhaeger et al. (2010) for WASP-18, r Boo, HD 179949, and v And, respectively; log(-^-) 
is calculated. 

a Estimated from our measured logLx < 27.5 erg s -1 and the fi HK — log (Lx/Lt, i) correlation presented in Mamajek & Hillenbrand 
(2008). See also Footnote 15. 

b The potential strength of magnetic star-planet interaction in the r Boo system is reduced due to the apparent tidal locking of the star 
to the planet; see §1. 

in HD 179949 to explain our non-detection for a simi- 
lar interaction efficiency, which seems implausible given 
their similar stellar spectral types. 

The lack of evidence for star-planet interaction in 
WASP-18 is unexpected given that we selected it as the 
system with the greatest predicted interaction strength 
(based on planetary mass and semi-major axis) and given 
the broad similarity in stellar properties between WASP- 
18 and other systems (notably HD 179949 and v And) 
for which star-planet interaction has been claimed. Star- 
planet interaction in WASP-18 would seem to be at best 
highly transient as it was not demonstrably present dur- 
ing our observations. For reference, Ca II H and K vari- 
ability has been described as phased with the planetary 
orbit in < 50% of the epochs at which HD 179949 and 
v And have been observed. Further, the timescale for 
a magnetic reconnection event in the simulations of Co- 
hen et al. (2011) is short, ~0.25 orbits. On the other 
hand, our spectroscopic coverage of WASP-18 extends 
over 8 complete planetary orbits, and the X-ray obser- 
vations cover many more (albeit at lower sensitivity per 
orbit). To the extent that star-planet interactions are 
transient, they should intuitively occur with greater fre- 
quency in extreme systems such as WASP-18, which con- 
trasts (if not definitively) with our findings. It might be 
alternatively suggested that, if our observations did not 
simply happen to take place within a period of relative 
quiescence, the efficiency of star-planet interaction for 
WASP-18 could be lower than has been identified in the 
past for less extreme systems. However, we emphasize 
again that WASP-18 is not notably distinctive in terms 
of stellar properties from HD 179949 and v And. (Below, 
we explore the possibility that the particularly massive 
and close-in planet in WASP-18 actually acts to suppress 
stellar activity.) 

Regardless of its underlying cause, this lack of ob- 
served variability in WASP-18 demonstrates that even 
extreme systems, arguably the best candidates to display 
planet-induced activity enhancements, challenge prevail- 
ing ideas concerning star-planet interaction. The ap- 
pealing prospect of calibrating star-planet interactions 
to estimate exoplanet magnetic field strengths would 
currently appear to require additional unambiguous evi- 
dence of such interactions occurring. 



3.2. Expected intrinsic activity in WASP-18 

We next consider the overall baseline level of chromo- 
sphcric and coronal activity in WASP-18, to provide con- 
text for the non-detection of planet-induced variability. 
WASP-18 shows atypically low Ca II H and K core emis- 
sion for its spectral class and inferred age (contrast with, 
e.g., HD 111456; Freire Ferrero et al. 2004). When stellar 
age is calculated from Ca II activity, main-sequence F, G, 
K, or M stars with inferred ages of 340-740 Myr (i.e., sim- 
ilar to the age of WASP-18, which is not based on Ca II 
activitjQ) are chromospherically "active" (following the 
definition of Henry et al. 1996), with characteristic values 
of -4.5 < R HK < -4.4 (from data in Wright et al. 2004). 
While stars with 0.48 < B - V < 0.52, as for WASP-18, 
display a wide range of activity, they are significantly less 
likely to have i? HK < -5.0 (19%, versus 37% for F, G, K, 
and M stars; from data in Wright et al. 2004). However, 
i? HK in WASP-18 is low; if estimated indirectly based on 
the L? HK — log (Lx/^boi) correlation presented in Mama- 
jek & Hillenbrand (2008), the measured X-ray limit for 
WASP-18 suggests i?n K ~ - 5 -°03 

The X-ray luminosity for WASP-18 is also lower than 
predicted. Nine nearby stars of spectral class F5V-F7V 
(i.e., matched to WASP-18) in a volume-limited ROSAT 
sample have X-ray luminosities of 27.5 < logLx < 
29.0 erg s"\ with mean and median values of log Lx — 
28.3 erg s _1 (from data in Schmitt & Liefke 2004). How- 
ever, WASP-18 is younger than typical main-sequence 
stars and should consequently be more X-ray luminous. 
While we are not aware of a specific comprehensive study 
of young F stars, in general the X-ray luminosity of 
F5+ stars modestly exceeds that of G stars (Schmitt 
1997) and Hyades G stars (hence of a similar age to 
WASP-18) have mean logL x ^ 28.8 erg s _1 (Stelzer 
& Neuhauser 2001; Preibisch & Feigelson 2005). Con- 
sequently, the XRT-derived limit of log Lx ~ 27.6 erg s _1 
marks WASP-18 as notably X-ray-faint relative to com- 
parable stars, with a log (Lx/^boi) ~ —6.2 similar to that 



14 The age of WASP-18 is estimated as 630±53° Myr by Hellier 
et al. (2009) based on stellar isochrones, and as 579_2so Myr by 
Brown et al. (2011) based on tidal interaction modeling. 

15 A recent high signal-to-noise Keck/HIRES spectrum of 
WASP-18, calibrated as described in Wright et al. (2004), gives 
_R HK = —5.15 (Howard Isaacson 2012, private communication). 



of the Sun. 

In summary, not only is there no apparent enhance- 
ment in chromospheric or coronal activity in WASP- 18 
that might be linked to interaction with the planet, but 
rather the above results, taken together, suggest that 
WASP-18 is unusually quiet for a young F6 star. It is 
possible that the stellar age is incorrect; the lcr error bars 
on the measurement given by Hellier et al. (2009) encom- 
pass > 1 Gyr, and the Brown et al. (2011) estimate de- 
pends upon complicated planet-star interactions. If the 
planet has acted to spin-up the star (recall the orbital 
plane and the stellar rotation axis are aligned), the cur- 
rent relatively rapid stellar rotation rate would not be 
reflective of true age. Stellar age has proven difficult to 
determine in some other hot-Jupiter systems; for exam- 
ple, Schroter et al.(2011) find that CoRoT-2a is X-ray 
bright and young but a late-K companion is X-ray un- 
detected, inconsistent with the inferred system age, and 
a similar situation may apply for HD 189733 (Pillitteri 
et al. 2011). If WASP-18 were 2-3 times older than cur- 
rent estimates, the intrinsic activity would still be some- 
what low but to a much less unusual degree. If the age 
is indeed correct, it is natural to consider whether the 
extremely massive and close-in planet could act to sup- 
press, rather than enhance, stellar activity. The tidal 
force exerted by WASP-18 upon its parent star is much 
greater than is typical even for hot Jupiter systems, and 
in fact Arras et al. (2012) note that the small but appar- 
ently non-zero eccentricity indicated by radial-velocity 
data is likely due to tidal fluid motion on the star. Mid- 
type F stars have shallow outer convective zones; perhaps 
the tidal pull is sufficient to repress dynamo activity in 
WASP-18 or comparably extreme systems (of which none 
are, however, currently known). This possibility could be 
assessed through magnetohydrodynamic simulations. 

3.3. An additional challenge to observations of 
star-planet interaction 

Here we explore the possibility that a hotspot rotat- 
ing upon the stellar surface can potentially mimic a 
shorter-period signature of star-planet interaction if only 
observed over a small fraction of the stellar rotation. 
Figure 5 shows simulated examples for HD 179949 and 
v And, which have stellar/planetary periods of 7/3.09 
and 14/4.618 days, respectively (Shkolnik et al. 2008). 
The half and full sinusoid models for planet-induced ac- 
tivity chosen for these examples are similar to those ap- 
plied by Shkolnik et al. (2005, 2008) to these systems, 
but here for an edge-on inclination and with arbitrary 
normalizations. (By construction the stellar hotspot am- 
plitude is taken as unity, and the relative amplitude of 
the interaction model is then greater for HD 179949 than 
for v And, also qualitatively similar to the published 
models.) The probability of the stellar and planetary 
phases co-aligning to mimic star-planet interaction for 
these particular models is not very large (e.g., ~3% of 
10000 random phase offsets for each example yield to- 
tal squared residuals less than 0.1), but with additional 
freedom to adjust the relative amplitudes, or to consider 
additional parameters (such as spot latitude or system 
inclination, within a limited range) or to choose from al- 
ternative model functional forms, or with a greater toler- 
ance for (apparent) outliers, this effect could potentially 
present a significant source of contamination within a 
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Figure 5. Examples of how a hotspot rotating with the 
star (blue line) could potentially mimic various models of 
star-planet interaction (black lines) over short observing runs 
(filled circles). The stellar hotspot phased with the planetary 
orbit is shown with open circles. See §3.3 for discussion. 

large sample of tested stars. For the particular cases of 
HD 179949 and v And it must be emphasized that several 
observing runs were conducted, and in some (but not all) 
of those runs variability of similar amplitude similarly 
phased with the planetary orbit was observed (Shkol- 
nik et al. 2005, 2008; Poppenhaeger et al. 2011a), which 
clearly decreases the probability of mistaking a stellar 
hotspot for one co-rotating with the planet. 

We note that, in general, searches for star-planet inter- 
action based on Ca II core variability can guard against 
this type of false positive at a given epoch through ob- 
taining coverage over at least one complete stellar ro- 
tation and over multiple planetary orbits (as was done 
with our coverage of WASP-18). If searches over mul- 
tiple epochs find some instances in which core variabil- 
ity phases with the stellar rotation and others in which 
it apparently matches the orbital period, a comparison 
of the amplitudes can also check whether the observed 
activity is likely to arise from distinct (intrinsic versus 
planet-induced) sources. It would be odd if the activ- 
ity in the "off" state (phased with stellar rotation rather 
than orbital period) had a similar amplitude to that seen 
in the "on" state, during which there is no obvious phys- 
ical reason that the star should otherwise go silent. It 
is not clear to us that the current data for HD 179949 
and v And can definitively pass this test, but see, e.g., 
Fares et al. (2012). In any case, the conclusion is that 
complete and extended coverage, preferably across mul- 
tiple epochs and multiple orbits per epoch, is essential 
for high-confidence detections of star-planet interaction. 

4. CONCLUSIONS 

We have carried out a sensitive search for planet- 
induced stellar activity within the extreme WASP-18 
system, selected as an ideal testbed for investigating 
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potential magnetic (or tidal) interactions between "hot 
Jupiters" and their host stars. High-resolution spec- 
troscopy of the Ca II H and K lines was conducted with 
the 6.5m Magellan Clay Telescope and contemporaneous 
X-ray monitoring was obtained with the Swift satellite. 
Our primary results are the following: 

1. The Ca II H and K cores do not show significant 
variability over ~8 d. Stacking residual spectra from 13 
visits into phase bins of 0.70-0.93, 0.99-1.19, and 1.21- 
1.41 does not show any significant structural changes 
that could be attributed to planetary influence. 

2. WASP-18 is not detected in a stacked 50 ks XRT 
exposure, constraining the X-ray luminosity to be unusu- 
ally low for a young F6 star, with logLx < 27.6. 

3. The lack of observed variability phased with the 
planetary orbit suggests any magnetic star-planet inter- 
action in WASP-18 must be transient, if present at all. 

4. The low level of chromospheric and coronal activ- 
ity is consistent with an intrinsically weak magnetic field 
on WASP-18, perhaps indicating that the inferred young 
stellar age is not reliable, or alternatively potentially re- 
lated to particularly strong planetary tidal effects. 

5. It is demonstrated that a stellar hotspot can po- 
tentially mimic star-planet interaction, for observations 
truncated to a short segment of the rotation period. 

6. Current ideas concerning star-planet interaction 
do not appear to be supported by the above results, 
therefore it may be optimistic at present to conceive of 
star-planet interaction as a robust estimator of exoplanet 
magnetic field strengths. 

Further high-quality Ca II H and K spectroscopy (as in 
Shkolnik et al. 2005, 2008) of previously monitored and 
newly discovered systems would refine understanding of 
the observational signatures of star-planet interaction. In 
addition, we are currently carrying out a Chandra sur- 
vey of solar analogs to check whether stars with close-in 
planets are systematically enhanced in X-ray luminosity, 
and this experiment has been designed to sidestep many 
of the selection biases that necessarily challenged previ- 
ous large-sample X-ray studies. That investigation will 
help further clarify whether the difficulties in establish- 
ing observational evidence of star-planet interaction in 
WASP-18 are anomalous or typical. 

We thank Mario Mateo for essential assistance with 
the Magellan observing run, Neil Gehrels and the Swift 
Science Team for approving and carrying out our Target 
of Opportunity observations of WASP-18, and an anony- 
mous referee for helpful comments. This paper includes 
data gathered with the 6.5m Magellan Telescopes located 
at Las Campanas Observatory, Chile. This research has 
made use of the Exoplanet Orbit Database and the Exo- 
planet Data Explorer at exoplanets.org. The Center for 
Exoplanets and Habitable Worlds is supported by the 
Pennsylvania State University, the Eberly College of Sci- 
ence, and the Pennsylvania Space Grant Consortium. 
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